function x = initialize_variables(N,M,V,C,min_range,max_range)
min = min_range;
max = max_range;
K = M + V + C;
k = V;
X = lhsdesign(N,k);
x = [];
range = [0 0 0 0 0 0 100 500;10 50 50 900 600 600 900 1200];
for o = 1
    x = [x round(range(1,o) + X(:,(o-1)*24+1 : o*24)*(range(2,o)-range(1,o)))];
end
for o=2:8
    x = [x round(range(1,o) + X(:,(o-1)*24+1 : o*24)*(range(2,o)-range(1,o)))];
end
for i = 1 : N
    x(i,V + 1: K) = evaluate_constraint(x(i,:),M,V);
end
    
end
